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We have carried out density-matrix-renormalization group (DMRG) calculations for the problem 
of one doped hole in a two-leg t — J ladder. Recent studies have concluded that exotic “Mott” 
physics — arising from the projection onto the space of no double-occupied sites — is manifest in 
this model system, leading to charge localization and a new mechanism for charge modulation. In 
contrast, we show that there is no localization and that the charge density modulation arises when 
the minimum in the quasiparticle dispersion moves away from tt. Although singular changes in the 
quasiparticle dispersion do occur as a function of model parameters, all the DMRG results can be 
qualitatively understood from a non-interacting “band-structure” perspective. 


A strongly correlated quantum system is one in which 
the interactions are at least comparable to the kinetic 
energy so weak-coupling, perturbative approaches can¬ 
not be justified. However, a key question is - under what 
circumstances does the behavior of such systems extrap¬ 
olate smoothly to the weakly interacting limit so that, 
at least at the phenomenological level, weak coupling in¬ 
tuitions can still be applied? There are certainly forms 
of broken symmetry, such as charge-density wave order 
in more than ID, which are at the very least unnatu¬ 
ral at weak coupling, and there can be still more exotic 
phases, especially those that support topological order 
and fractionalization, which have no weak-coupling ana¬ 
logues. What about the important case of a doped Mott 
insulator? It has been argued by many authors that there 
is an additional quantity, sometimes referred to as “Mot- 
tness”, which through the effect of the constraint of no 
double-occupancy produced by a strong local “Hubbard 
7/,” can invalidate the quasiparticle picture and preclude 
the adiabatic continuation to the weakly interacting ref¬ 
erence state that underlies Fermi liquid theory. 

The idea that the quasiparticle picture fails qualita¬ 
tively hasgained strong support from a set of papers by 
Zhu et in which extensive numerical experiments 

have been carried out using the density matrix renormal¬ 
ization group (DMRG)I^ on a set of t — J ladders. It has 
long been thought that the undoped two-leg t — J ladder 
is adiabatically related to a band insulator, and a num¬ 
ber of early exact diagonalizatioiP and quantum Monte 
CarlcP studies supported the idea that doped holes form 
conventional quasiparticles. In striking contrast, Zhu et 
al reported that a doped hole in a two leg ladder local¬ 
izes at large length scales, a finding that is incompatible 
with Bloch’s theorem for any quasiparticle state. Similar 
localization was reported on three and four leg systems, 
although the data is less extensive. Zhu et al proposed an 
explanation for this behavior based on considerations of 
“hole phase-strings” and a new type of “Weng statistics.” 
It has been further proposed,!^ that this new paradigm 
can account for a wide range of phenomena in doped Mott 
insulators, including stripe formation in the cuprates. 

In this paper, we have focussed on the two-leg t — J 


ladder with one doped hole. We have carried out DMRG 
calculations to extract the ground-state properties of lad¬ 
ders of length up to L = 1000, and time-dependent 
DMRGP^^ (tDMRG) calculations on ladders up to L = 
120 to obtain unprecedentedly complete information con¬ 
cerning the dynamical one-hole Green function, G. Fol¬ 
lowing Zhu et al we have considered a range of values 
of the parameter a, the ratio of the hopping matrix el¬ 
ements and the exchange couplings on the legs and the 
rungs of the ladder. In contrast to them, we find that 
the one hole state is never localized. On the other hand, 
we corroborate their discovery that a notable change in 
the character of the one-hole state occurs at a critical 
value of a = Q^c ~ 0.68; in particular the quasiparticle 
effective mass diverges as a ^ o^c. However, this sin¬ 
gular behavior does not imply the existence of a phase 
transition, as changes in the properties of a single doped 
hole do not reflect changes in the thermodynamic state of 
the system. Indeed, we show directly from the structure 
of G that the quasiparticle is well defined for a on both 
sides of Q^c, that there is no “spin-charge separation,” and 
that the quasiparticle weight, Z{a), is always substan¬ 
tial. Indeed, all the properties of the low energy one hole 
states can be adiabatically related to those of a single 
hole in a non-interacting “band” insulator - the singular 
changes reflect a shift of the ground-state sector from a 
Bloch wave vector k = tt for a < ac to k = ko{a) < tt 
for a > ac. The divergent effective mass dramatically 
reflects a point at which the minimum of the quasihole 
dispersion, e{k), shifts away from tt. 

In this paper we will study the 2-leg t — J — a model 

H = — ~ (1) 

{hj),cr {i,j) 

Here (ij) indicates nearest-neighbor sites with tij = t 
and Jij = J on the rungs, and Uj = at and Jij = aJ on 
the legs, creates an electron on site j with spin po¬ 
larization cr, the spin operator on site j is Sj, the charge 
is Uj = action of the Hamiltonian 

is restricted to the Hilbert space with no doubly occu¬ 
pied sites, Uj = 0, I. The index i = with ly = I 
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and 2 denoting the two legs and lx runs from 1 to L. 
This is the same realization of the t — J model that was 
studied by Zhu et al. for a range of a with J/t = 1/3. 
They gave a quasiparticle interpretation to their results 
for Of < ac ~ 0.7, but they identified a transition at 
a = ac^ such that, among other anomalies, for a > ac 
and ladders of length L > 100, they reported localization 
of the charge in a region of width ^ ^ 100 < L. 

Our ground state DMRG calculations were fairly stan¬ 
dard, the main exception being that an unusually large 
number of sweeps were needed for the one hole ground 
states. All the calculations reported here were performed 
using the ITensor library (http://itensor.org). A suffi¬ 
cient number of states, roughly 200-400 for the one hole 
case, were kept to limit the truncation error per step 
to ~ 10“^^. For each system, first the ground state for 
the undoped system was obtained, with four sweeps giv¬ 
ing high-accuracy convergence, and this matrix product 
state 10) was stored. We then applied the operator 
where jo is a site at the center of the system, creating 
a one hole state with the hole localized in the center. 
Sweeps were then carried out, resulting in a set of ever 
better approximate one-hole groundstates, |'0(s)), where 
s indicates the number of sweeps. At each sweep we made 
diagonal measurements of the energy and the density on 
each site, as well as off-diagonal measurements of the hole 
amplitude, F{j,s) = ((/)|c]^|V»(s)). 

Figure 1(a) shows the spreading of the density in a 
1000 X 2 system versus sweep with a = 1. Here the hole 
density for site j is (j) = 1 — rij ; the figure shows the 
rung hole density nh{lx) = The density 

continues to spread out as the sweeps progress. (Note, 
to facilitate comparisons with previous results, we have 
eschewed tricks that could be used to accelerate conver¬ 
gence to the true ground state, such as starting with 
a delocalized hole as the initial state.) The inset in 
Fig.j^a) shows the full width at half-maximum (FWHM) 
of the charge density profile for a = 1 ladders of different 
lengths L. This value of a is greater than ac and places 
the system in the region where Zhu et al reported lo¬ 
calization. However, as seen in the inset, we find that 
the FWHM scales as L. The saturation of the FMHW 
reported by Zhu et al. in Fig. 2c of Ref [4] appears to be 
an artifact of their calculation which arises from limiting 
the number of DMRG sweeps. In fact, as shown in Fig. 
4c of [4], they, too, find the charge density extends over 
a 200x2 ladder when the sweep number is increased. 

Figure 1(b) shows a correlation function 
{S^{lx^ly)Tih{jo)) which measures the spin profile 
when a dynamic hole is on site jo; here jo = (200, 2) on a 
400 X 2 ladder. With this correlation function shown on 
a log scale as a function of distance lx along the ladder, 
the exponential confinement of the spin and charge is 
apparent in the linear lx dependence. A linear fit gives 
a decay length of ^ = 3.14 for a = 1; this matches 
closely with previous results of 3.19(1) for the spin-spin 
correlation length in the undoped ladder(In contrast, 
Zhu et al. reported that a similar correlation function 


decayed as a power law ioi a > acj 




FIG. 1: (a) The density on each site for o = 1 on the cen¬ 
tral portion of a 1000 x 2 system for the indicated number of 
sweeps. The inset shows the full width at half maximum of 
the density on a set of smaller lattices which were converged 
in the number of sweeps, (b) A correlation function which 
measures spin-charge correlations, showing that the spin de¬ 
grees of freedom are exponentially localized close to a dynamic 
hole, for a = 0.5 and a = 1. For a = 1, the red line shows a 
linear fit to the data. 

To obtain the one particle spectral function, instead 
of evolving \jj) with DMRG sweeps, we evolve it in real 
time, obtaining the state 

IV’(^)) = (2) 

After each time-step, the Green function, 

= (3) 

(defined here without the usual i prefactor) was measured 
for all sites j, where Eo is the ground-state energy of the 
undoped ladder. As time evolves, the wavepacket spreads 
out. We always stop the simulation at a time t^ax before 
the packet reaches the edges of the system. Thus any 
finite size effects are completely negligible, arising only 
from the undoped state, which has a correlation length 
that is very small compared to L. Other sources of error 
are the finite tmax^ finite truncation error, and finite size 
of the time steps. Using time steps r = 0.05 — 0.1, we 
found the time step error was small enough to have no 
visible effects on any of the figures below. To measure 
and control the other two errors, we varied the number 
of states kept (up to m = 2000) and the maximum time 
(up to tmax = 100). Any errors in the results we show 
primarily appear as slight broadenings of the spectra, and 
have no impact on our conclusions. 

The ladder is symmetric under reflection symmetry 
which interchanges the two legs; correspondingly, the 
one-hole states can be classified by their symmetry, A = 
±1, under reflection. Similarly, the Bloch wave-number is 
a good quantum number. Thus, to interpret the results 
physically, we perform the Fourier transform of G(j, t) 
with respect to time (using G{j,—t) = G{j^tY) and 
position along the ladder, projected onto the space of 
states of agiven reflection symmetry using both linear 
predictiorP^ and windowing to deal with a finite tmax- 
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The real part of this quantity is the spectral function 
A{k^uj)^ which is shown for A = +1 in Fig. [^a) for the 
case a = 1. The Supplementary Information section con¬ 
tains a further discussion of the tDMRG and figures of 
A{k^uj) for more values of a. 

The spectral weight is characterized by a sharply de¬ 
fined dispersing pole separated by a gap of order J 
from a quasi-particle-magnon continuum. For a = 1, 
the minimum in the quasi-particle dispersion occurs at 
kmin ~ 2.01 = 0.6407r. A slice of the spectral weight for 
a = 0.7 (just above ac ~ 0.68) at kmin ~ 2.85 = O.OOTtt 
is plotted versus uj in Fig. 2(b). The dispersion of the 
pole in the quasi-particle spectrum versus k for several 
values of a is shown in Fig. 2(c). As a increases beyond 
5 kmin moves away from tt and at large values of a 
approaches 7r/2. 
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FIG. 2: (a) Spectral weight function A(k^uj) in the ground- 
state (A = H-l) reflection parity sector for the t-J ladder with 
0 = 1 , obtained with tDMRG color indicating the value of 
A(k,(ju). We work in energy units where t — 1. (Results 
for odd reflection parity, A = — 1, are shown in the Supple¬ 
mental Section.) (b) A{k,uj) near the quasiparticle peak for 
a = 0.7 at ko/iT = 0.907. The gap to the start of the con¬ 
tinuum spectrum is of order J. (c) Quasiparticle dispersions 
for a = 0.5, 0.7, 0.9, obtained from tDMRG. The stars show 
the values of ko and sq obtained from separate ground state 
DMRG calculations. 

For a given value of a, the minimum hole energy sq 
and the corresponding wave vector ko can be determined 
from the dispersion of the peak in A{k, uj). Alternatively, 
for a given value of a, the energy sq and wave vector ko 
can be determined directly from our ground state DMRG 
calculations. The energy minimum so for a given value of 
a is equal to the difference in the one hole and zero hole 
ground state energies. The wave vector ko associated 
with the one-hole ground state can be determined from 
the peak in the spatial Fourier transform of F(jf, s), which 
sharpens as the sweep number s increases. Plots of So and 
ko versus a are shown in Fig. 3. 



FIG. 3: The one quasi-hole properties as a function of a: The 
figure shows (a) Z and £o, and (b) ko and m*. These results 
were obtained from 400 x 2 systems from measurements as 
the hole spread out with successive sweeps. 

Similarly, while m* can be extracted from the curva¬ 
ture of the quasi-particle dispersion around ko and Z 
can be obtained from a frequency integration of A(/co, cj), 
both of these quantities can be directly determined with 
higher accuracy from the ground state DMRG calcula¬ 
tions. An estimate of the quasi-particle spectral weight 
Z is given by 

Zis)=J2\Fij,s)f (4) 

j 

We find that this estimate converges very rapidly with 
the number of sweeps s, much more rapidly than the hole 
spreads out. As the DMRG sweeps continue, the energy 
E of \'ip{s)) converges towards that of the one-hole ground 
state with a correction that varies as (8m*(x^))“^. Here, 
is the variance of the position of the hole, determined 
from {nh{lx))‘ By plotting E versus with each 

point corresponding to a different sweep, one can obtain 
an estimate of m*. In addition, one can increase the 
accuracy of the estimate for Z for the infinite ladder by 
extrapolating Z versus For a = 1, for example, 

we obtain Z = 0.34067(1). Plots of Z and m* are shown 
in Fig. 3. As seen in this figure, there is a sharp change 
in the quasi-particle character that occurs at ac = 0.68. 
There are kinks in the slopes of So^ko and Z and the 
curvature of the quasi-particle dispersion vanishes giving 
rise to a divergence in the effective mass. The shift in ko 
away from tt gives rise to the oscillations in the charge 
density, as has been previously noted by Zhu et al, which 
are found to occur at wave-number 2ko. 

Since the one hole state has a well defined quasi¬ 
particle spectral weight, many properties that are mea¬ 
surable in numerical experiments on systems with large 
but finite L can be understood in terms of the simpler 
problem of one-hole on a 2-leg band insulator. Gentral 
to this understanding is the quasi-particle dispersion re¬ 
lation which determines the values of /c = ±/co at which 
e{k) is minimized, and the dependence of the hole-energy 
near this point, e{k) = Eo E eo E {k — ko)‘^/2rrE + ..., 
where rrE is the effective mass. In order to minimize its 
zero-point energy on a ladder of large but finite length 
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L, the one-quasiparticle ground-state will always spread 
to fill the extent of the ladder, 

^ sin(7rn/I/) cos(A:on — 6>), (5) 

where 0 = kQL/2. The minimum in the ground state 
energy of the one hole state is 

s{L) = Eq +50+ TT^/ (2m ^+ .... (6) 

Since integrating out the gapped spin degrees of freedom 
inevitably renormalizes the bare dispersion, for compar¬ 
ison purposes we consider a non-interacting model with 
band structure 

E{k) = —— 2t|| cos{k) — 2tj| cos(2/c) (7) 

in which A = ±1 correspond to the valence and con¬ 
duction bands, respectively, all the t’s are assumed non¬ 
negative and the rung hopping parameter t± to be suf¬ 
ficiently large compared to the near-neighbor and next- 
near-neighbor leg hopping parameters ty and tj| that the 
undoped system has an insulating gap. This dispersion 
is similar to that shown in Fig. 2c. The parameter that 
plays a role analogous to a is d = 4tj|/t||; for 0 < d < 1, 
the top of the valence band occurs at /c = tt, while for 
d > 1, the top of the valence band occurs at k = ±ko 
where cos(/co) = — 1/d. The critical dependences of 
So = —E{ko) , ko and m* on d can be readily derived 
from the band dispersion Eq. (7). 

[sq Et±] _ f -(4-d)/2 for d<l . . 

t|| \ —(2 + d?)l2a for d > 1 ^ ^ 


_ J 1/2 for d < 1 . . 

t|| da ~ \ (2 - 6?) 126? for d > 1 ^ ' 

{ 0 for d < 1 

A/2(d — l)/d for 1 (d — 1) > 0 (10) 

7r/2 — 1/d for d^l 


and 

*_J_r fora<l , . 

2i|| \ — 1) ^ for a > 1 ^ ' 

The qualitative features observed in the evolution of the 
one-hole state of the t — J — a model as a function of 
a are reflected in the band model as a function of d. 
i) The one-hole energy sq has a non-analytic change in 
slope at d = dc given by Eq. [9]. ii) The vector ko{a) 
has a square-root singularity at d = dc as given by Eq. 
[10], and 2ko determines the oscillations of the charge 
density, iii) The effective mass m*(d) diverges linearly 
upon approaching dc from both sides as given in Eq. [11]. 

In the Supplemental Information we make this con¬ 
nection formal: We define at — J-Hubbard model Hamil¬ 
tonian that in one limit is equivalent to the t — J — a 
model of Eq. 1, and in another limit represents a non¬ 
interacting band-insulator, with the band structure given 
in Eq. (7). Our DMRG results establish that there is 
no gap-closing and so no barrier to adiabatic continuity 
upon reducing the model to one of decoupled rungs in 
the a = 0 limit. In this limit the interactions can be 
adiabatically set to 0, again without any gap closures. 
Einally, in the solvable non-interacting limit, we restore 
the hopping matrix elements along the ladder, ty and tj|, 
still without encountering any gap closures. (The final 
two steps are readily studied analytically.) This analysis 
constitutes a proof that the low energy one-hole states of 
the t — J—a model are adiabatically connected to those of 
a non-interacting band insulator which holds regardless 
of the value of a in the entire range we have studied. 
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Appendix A: Supplemental material 


has a non-interacting limit, 


The fact that the low energy one-hole states of the two- 
leg t — J ladder can be adiabatically connected to the 
corresponding states of a two-leg non-interacting semi¬ 
conductor is established by explicit construction. In ad¬ 
dition, details of the time dependent DMRG calculations 
are given, and additional plots of the spectral function 
are presented. 


1. Adiabatic continuity to the noninteracting limit 


The t — J Hubbard model: Consider the Hamiltonian 
for electrons in a two leg ladder with sites labeled by the 
leg index r = d and rung index j: 


~ y ] ^\\^],T,a^j+‘2,T,a + H.C. 

j,ra 

~ + H.C. 

ho- 


+J± ^ 


J,r 




^j,u ’ ^j,d 




(Al) 


J,o 


where ty and tj| are the first and second neighbor hopping 
along the ladder, t±_ is the hopping between rungs, Jy 
and J± are the corresponding exchange couplings, cj ^ 
creates an electron with spin-polarization a on site (j, r), 


~ T 




a,a' 


(A 2 ) 


is the spin and 


^ 3,0 — ^J,r,cr^i,T,cr (^3) 

a 

is the electron density on site In contrast to the 

t — J model, this hamiltonian acts on the full fermionic 
Hilbert space in which there is no constraint on double¬ 
occupancy sites, although this constraint can be obtained 
dynamically by taking the limit in which the on-site Hub¬ 
bard repulsion U tends to oc. Thus, the t — J model is 
simply the U ^ oo limit of this model; specifically, the 
version of the t — J model studied by Zhu et al (Eq. ( 1 ) 
of our paper) is 

Ht-j = lim 0, t, (aJ, J, I/). (A4) 

U^oo 

with t = ty, J = Jy, and a = = Jy/Ja- On the 

other hand, in contrast to the t — J model, this model 


Hnonit,t',t±) = 0,0,0). (A5) 

Non-interaeting limit: The band-structure of the non¬ 
interacting two-leg ladder described by Hnom is trivially 
obtained using Bloch’s theorem and taking advantage of 
the reflection symmetry which exchanges the two legs. 
The dispersion relation for this problem is 

e^{k) = —ytx — 2 ty cos{k) — 4tj| cos^(/c) - 1 - 2 tj| (A 6 ) 

where 7 = ±1 is the reflection symmetry and k is the 
Bloch wave-number. We impose the condition the system 
be insulating when there is one electron per site by con¬ 
sidering \t± \ is sufficiently large compared to |ty | and |tj| |, 
so that there is a gap in this spectrum. To be explicit, 
we further restrict consideration to the case in which all 
the t’s are non-negative. The parameter that plays a role 
analogous to a is d = tj|/4ty; for 0 < d < 1, the top 
of the valence band occurs at /c = tt, while for d > 1, 
the top of the valence band occurs at /c = ±ko where 
cos(/co) = 1/d. 

Adiabatie eontinuity: In a quantum system with a gap, 
the notion of adiabatic continuity can be given a precise 
definition - the states of two systems are adiabatically 
connected if it is possible to continuously deform the 
Hamiltonian in such a way that the gap never closes in 
turning it from that of the initial to the final system. If 
we restrict ourselves (as we often do) to adiabatic paths 
that preserve certain symmetries, then two states with 
different symmetry related quantum numbers can never 
be adiabatically connected. Conversely, we can study 
adiabatic continuity within a given subspace of Hilbert 
space specified by these quantum numbers, even if some¬ 
where along the path there might occur a region where 
the absolute ground-state lies in a different subspace.® 
An adiabatie route from Hf-j to Hnon- Here we vary 
the parameters in , tj|, Ty, t/) to trace an adi¬ 

abatic path from Ht_j to H^on always preserving trans¬ 
lational (with periodic boundary conditions), reflection, 
spin rotational, and gauge (number conservation) sym¬ 
metries: 

Ht-j = H{at^ OH] = oc) (A7) 

—^ 77(0, 0 , tx; 0 , 1/ = oc) 77(0, 0, t±; 0, 0; 0) 

^ H{tii,-t\^,t±-,0,0;U = 0 ) = Hnon{t\\,t\^,t±). 

The single-arrow represents steps for which the existence 
of a non-zero gap along the entire path and hence the pos¬ 
sibility of adiabatic evolution has been established using 
DMRG results for the t — J ladder, while the double¬ 
arrows represent steps that can be justified analytically. 

i) In the first step, the system is deformed into a set 
of decoupled rungs. That there is a (spin) gap in the ex¬ 
citation spectrum of the undoped ladder over the entire 
pertinent range of a is well established by the present 
DMRG calculations as well as by those of Zhu et aW^. 
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For one doped hole, we first restrict attention to the sub¬ 
space with k = ko and X = 1 - i.e. the sector which 
contains the one-hole ground-state for the particular ini¬ 
tial value oi a = ao being considered. (For o^o < olc this 
is /co = TT, whereas for > ac this is an appropriate 
smaller value of k^.) What is apparent from our DMRG 
results is that for the entire range a we have studied, 
there is a gap of order J separating the quasi-hole state 
from the multi particle continuum. Notice that although 
in the case that < ac^ the one-hole ground-state re¬ 
mains at /c = TT for the entire range of a between a = 
and (a = 0, for o^o > olc^ the ground-state sector changes 
as a varies from a = to a = ac- This does not, how¬ 
ever, act as a barrier to adiabatic evolution, since as long 
as we maintain translational symmetry, we are free to re¬ 
strict our attention to the subspace with k = ko{ao) for 
the entire process. 

ii) Once the system consists of decoupled dimers, the 
spectrum can be readily computed analytically and it 
is easy to see that it is possible to simultaneously de¬ 
crease the values of U and J_l, without ever closing the 
gap, to the point at which the system consists of non¬ 
interacting electrons confined to the bonding states on 
each rung. Note that throughout this portion of the evo¬ 
lution, the one-hole ground-state is 2L fold degenerate, 
but in any sector specified by Bloch wave-vector, /c, and 
spin polarization, <j, the one-hole ground-state remains 
non-degenerate. 

iii) Once all interactions have been quenched, it is sim¬ 
ple to compute the band-structure for arbitrary t_L, ty, 
and tj|. For sufficiently large the semiconducting 
gap of the undoped system the gap at fixed k in the 
presence of one doped hole remain non-zero throughout 
this process. If in the starting Hamiltonian, a < ac, 
then we can insure that the one-hole ground-state oc¬ 
curs in the appropriate ko = tt sector by ending with a 
value of |tj|| < |t|||/4. If the starting Hamiltonian has 
a > ac, then by ending the adiabatic evolution with 
tj| = —1||/4 cos(/co), we reach a situation in which the 
lowest energy one-hole state occurs at k = kg (a), i.e. in 
the same sector of Hilbert space as in the initial Hamil¬ 
tonian. 

This constitutes the proof that the two-leg t — J ladder 
is adiabatically connected to a band-insulator, both for the 
undoped system and in the presence of one doped hole. 


2. Time dependent DMRG and spectral functions 

The tDMRG results were obtained using a Trotter 
decompositiorPM^, applying only nearest neighbor gates, 
using a reordering of the sweep path through the lattice 
to make this possible. The initial path was chosen as 
{x,y) = ( 1 , 1 ), ( 1 , 2 ), ( 2 , 2 ), ( 2 , 1 ), (3,1), .... Along this 
path, only bonds on leg 1 are not nearest neighbor. We 
alternate this path with (1, 2), (1,1), (2,1), (2, 2), (3, 2), 
..., where leg 2 bonds are not nearest neighbor. Write 


the Hamiltonian as H = Hi -\- H2 where 
Hi = —-Hrungs T -f^leg—1 

and 

H 2 = “-^rungs 4“ -^leg—2 

We perform a half sweep, applying bond time evolution 
operators exp(—irH'bond), applied only on the terms in 
Hi, where r is the time step. Then, we perform a half 
sweep which switches the path to the second ordering, 
consisting of applying swap operators on each run^. Af¬ 
ter this, we perform a half sweep using the terms of H 2 . 
Then this entire sequence of three half sweeps is applied 
entirely in reverse, returning the sites to their original 
order. The reversal also cancels the lowest order Trotter 
error, resulting in an overall Trotter error of O(t^) (per 
unit time), in a time-step that progresses by 2 r. 
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FIG. 4: Spectral function for a = 0.5 for A = -hi. The 
maximum number of states was m = 500. 

A Fourier time-space Fourier transform of G{j, t) yields 
A{k,uj), using both linear predictiorP^ and windowing to 
deal with a finite tmax- As time evolves, the wavepacket 
spreads out. We always stop the simulation before the 
packet reaches the edges of the system. Thus any fi¬ 
nite size effects are completely negligible, arising only 
from the undoped state, which has a very short correla¬ 
tion length. Other sources of error remain, namely finite 
maximum time, finite truncation error, and finite time 











step. Using r = 0.05 — 0.1, we found the time step error 
was small enough to have no visible effects on any of the 
figures below. 
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FIG. 6: Spectral function for a = 0.7 for A = +1. The 
maximum number of states was m = 500. 


FIG. 5: Spectral function for o = 0.5 for A = — 1, for the 
same run as in the previous figure. 

As t increases, the entanglement of the state increases. 
If we keep a variable number of states, specifying a par¬ 
ticular truncation error at each step, then the number of 
states will increase as time increases. For example, with 
(a = 1 if we specify a truncation error of 10“^, the num¬ 
ber of states kept rises as a rapidly increasing function 
that reaches m = 3000 at about t ^ 14 — 15. (The en¬ 
tanglement growth is smaller for smaller a.) There are 
several ways to deal with this entanglement increase: we 
discuss three approaches. 

1) One can stop the simulation when m reaches a cut¬ 
off, e.g. stopping at tmax ^15 when m = 3000 for a = 1. 
One can rely on the linear prediction to extend tmax be¬ 
fore Fourier transforming. We did not follow this ap¬ 
proach. 

2) One can specify a maximum m to overrule the spec¬ 
ified truncation error, giving a larger truncation error for 
larger times. An advantage of this method is that the 
increased truncation somewhat resembles a windowing 
function, in that it reduces G{t) as t increases. A win¬ 
dowing function must be applied anyway, so this is not 
a very serious error. The decrease in G{t) is not uniform 
across frequencies-the higher energy states are more en¬ 
tangled, and their amplitude decreases more rapidly than 
the low energy states. This allows good resolution of the 
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FIG. 7: Spectral function for o = 0.7 for A = —1, for the 
same run as in the previous figure. 

























quasiparticle part of the spectrum, and we have generally 
followed this. (One must not “fix” the normalization of 
the wavefunction after the truncation error-this increases 
the amplitudes of the low energy part of the spectrum to 
fix the loss at high energies, producing poor results.) One 
can vary m, tmax, etc. and check for convergence of the 
results. We have mostly followed this procedure. 
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FIG. 8: Spectral function for ol — 1.0 for A = +1, keeping 
m — 2000 states. 

3) One can evolve in imaginary time a fixed distance p 
(not an actual temperature), say /3 ^ 1, before starting 
the real time evolution. This diminishes the high energy 
parts of the state, and the low energy part that is left 
has lower entanglement growth. After the simulation, 
and after the standard linear prediction, windowing, and 
Fourier transforming, one corrects for the initial imagi¬ 
nary time evolution by replacing A{k^uj) by A{k,uj)e^^. 
This method works quite well. For large cj, the results 
can be poor, because errors are amplified, but it gives a 
well controlled way to zoom in on the low energy part 


of the spectrum with high accuracy. This approach can 
be combined with method 2), using a maximum m. This 
method was used to obtain the spectrum shown in Fig. 
2(b). 

We present here the spectral functions for several val¬ 
ues of a, all obtained with method 2). The low energy 
quasiparticle bands are very reliable, and their finite fre¬ 
quency width is a consequence of finite tmax- The higher 
energy parts have broad features and also smaller am¬ 
plitude features, e.g. subtle color variations. The broad 
features and distribution of spectral weight are very reli¬ 
able, but it can be hard to tell if some of the small ampli¬ 
tude high frequency features are artifacts due to noise or 
ringing, without further study, comparing spectra with 
different accuracy parameters, which we have not done 
very thoroughly. 
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FIG. 9: Spectral function for o = 1.0 for A = — 1, for the 
same run as in the previous figure. 


^ The thermodynamic definition of adiabaticity is somewhat 
more complicated, involving the possibility of a quasi-static 
evolution at constant entropy. Since the one-hole state is 
not a sensible notion in the thermodynamic limit, and as in 
any case the quantum mechanical definition of adiabaticity 
is more restrictive than the thermodynamic one, we deal 


here solely with the former. 
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